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ABSTRACT 

We introduce a stellar velocity survey of dwarf spheroidal galaxies, undertaken using the Michi- 
gan/MIKE Fiber System (MMFS) at the Magellan/Clay 6.5 m telescope at Las Campanas Observa- 
tory. As of 2006 November we have used MMFS to collect 6415 high-resolution (R = 20000 - 25000) 
spectra from 5180 stars in four dwarf spheroidal galaxies: Carina, Fornax, Sculptor and Sextans. 
Spectra sample the range 5140 — 5180 A, which includes the prominent magnesium triplet absorption 
feature. We measure radial velocity (RV) to a median precision of 2.0 km s _1 for stars as faint as 
V ~ 20.5. From the spectra we also are able to measure the strength of iron and magnesium ab- 
sorption features using spectral indices that correlate with effective temperature, surface gravity and 
chemical abundance. Measurement of line strength allows us to identify interloping foreground stars 
independently of velocity, and to examine the metallicity distribution among dSph members. Here we 
present detailed descriptions of MMFS, our target selection and spectroscopic observations, the data 
reduction procedure, and error analysis. We compare our RV results to previously published mea- 
surements for individual stars. In some cases we find evidence for a mild, velocity-dependent offset 
between the RVs we measure using the magnesium triplet and previously published RV measurements 
derived from the infrared calcium triplet. In companion papers we will present the complete data sets 
and kinematic analyses of these new observations. 

Subject headings: galaxies: dwarf — galaxies: kinematics and dynamics — (galaxies:) Local Group 
— techniques: radial velocities 



1. INTRODUCTION 

The dwarf spheroidal (dSph) satellite companions of 
the Milky Way (MW) contain a wealth of information 
regarding the formation and ongoing evolution of the Lo- 
cal Group. The ~ 15 known MW dSphs are the nearest, 
smallest, and faintest known galaxies ( 0.5 < R/kpc < 1; 
10 4 < L/L0 < 10 7 ; 20 < D/kpc <; 250; lMateolll998l and 
refere nces therein: IWillman et alJl2005t iBelokurov et al.l 
2006). From their resolved stellar populations one can 
study in detail dSph morphologies, star formation histo- 
ries, chemical abundances, and kinematics. Given their 
surface brightness profiles, the stellar velocity dispersions 
(<r v ~ 10 km s _1 ) of these galaxies indicate a larger mass 
than could be contributed by any reasonable equilibrium 
stellar population. Applying the virial theorem, one de- 
rives for the MW dSphs masses of ~ 10 7 — 10 s M© and 
mass-to-li ght (M/L) r atios ranging from five to several 
hundred (Mateo 199a and references therein). To the 
extent that the virial theorem is applic able despite the 
tidal field exerted by the Milky Way (|Piatek fc Prvorl 
Il995t lOh efail 119951 ). the derived M/L ratios place at 
least some dSphs among the most dark-matter domi- 
nated galaxies known, generally with negligible baryonic 
contribution to their total mass. In this sense the MW 
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dwarfs offer an ideal testbed with which to constrain the 
behavior of dark matter at perhaps the most difficult 
scale to model — the lower extremum of the halo mass 
function. 

Efforts to characterize the internal kinematics of dSphs 
have expanded in s cope over the pas t quarter-century. In 
a pioneering study, lAaronsonl (|1983h considered the pos- 
sible dark matter (DM) content of dSphs based on the 
radial velocity (RV) dispersion of just three Draco stars. 
The need to invoke DM to describe dSph kinematics was 
placed on firmer footing afte r samples of t ens o f stars 
ilAaronson &~Q lszcwski 119871: IMateo et al.lli"99ll 119931: 
i Suntzeff et al.lll993t lOlszewski et al.l 1 19951: Queloz et al.l 
19951; IHargreaves et all 11994 Il996bi IVogt et al.l 119951 : 
IMateo et al.l [1998) confirmed Aaronson's result. RV 
samples for nearly 100 stars ac c umulated over several 
epochs (lArmandroff et al I 119951 : lOlszewski et al.l 119961 : 
IHargreaves et al.lll996al lbD showed that the magnitudes 
of dSph velocity dispersions were not unduly influenced 
by the presence of binary stars or motion within the stel- 
lar atmospheres. Recent RV samples for more th a n 100 
stars per dSph dKlevna et alj [200ll 120021 l200l l200l 
Tolstoy et alt [2 004; M unoz et all 120061 ; iBattagli a et all 



2005 iWalker et al.l 2006a b: Koc h et al.j|2006h . yield the 



general result that dSph velocity profiles appear isother- 
mal over most of the visible faces of dSphs, though there 
is some dispute o ver the velocity behavior near the nom - 
inal tidal radius (|Klevna et al.ll2004t iMunoz et al.ll2005h . 
Moreover, the growing number of measurements for stars 
at and even well b e yond the nomina l tidal radius (e.g., 
IMunoz et all 120061 ; [Sohn et alj 120061 ; Mateo et al. in 
preparation) has re-energized debate about the degree 
to which external tidal forces influence dSph stellar kine- 
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matics. 

Studi es bv iMerritt et alJ (|1997| ): IWilkinson et all 
(|2002l) ; IWang et all (|2005l ) have shown theoretically 
that velocity samples for ~ 1000 stars are capable 
of distinguishing among spherical mass models while 
relaxing some of the assumptions (e.g., constant mass- 
to-light ratio and velocity isotropy) characteristic of 
early analyses. The present observational challenge 
thus is to acquire data sets consisting of > 1000 stel- 
lar RVs with spatial sampling that extends over the 
entire face of the dSph. In order to resolve narrow 
dSph velocity distributions one must work at high 
resolution. Stellar targets sufficiently bright (V< 21) 
for measuring radial velocity and chemical abundances 
are generally confined to the dSph red giant branch. 
The development of multi-object fiber and/or slit- 
mask spectrographs with echelle resolution has fueled 
a recent boom in this field. Recent and ongoing 
work u sing VLT /FLA M ES/G IRAFFE (ITolstov et all 
120041: IBattaglia et~aT] 120061: IMunoz et all I2006D . 
WHT/WYF FOS (iKlevna et all 120041) . Keck/DEIMOS 
(jSohn et al.ll200l iKoch et alj|2006j) . MMT/Hectochelle 
(Mateo et al. in preparation) has yielded RV samples for 
hundreds of stars in several MW dSphs. These studies 
have produced several result s that challenge previous 
assumptions. ITolstov et al.l |2004) find evidence for 
multiple stellar populations in Sculptor that follow 
distinct distributions in position, metallicity, and RV. 
IBattaglia et alJ (|2006l) produce a sim ilar result in a 
study of Fornax. IMunoz et al.l (|2006l ) discover several 
widely separated Carina members, which they interpret 
as evidence of tidal streaming; they also attribute a 
secondary clump in Carina's velocity distribution to 
contamination by stars once associated with the Large 
Magellanic Cloud. Mateo et al. (in preparation) and 
Sohn et al. (2006) detect a kinematic signature of tidal 
streaming from the outer regions of Leo I. It is becoming 
clear that dSphs are complex systems that can be 
characterized adequately only with data sets of the sort 
coming from large surveys. 

Here we introduce an independent dSph kinematic sur- 
vey undertaken with a new instrument, the Michigan- 
MIKE Fiber System (MMFS). MMFS provides multi- 
object observing capability using the dual-channel Mag- 
ellan Inamori Kyocera Echelle (MIKE) spectrograph at 
the Magellan 6.5m Clay Telescope. Ours is the only 
large-scale dSph survey to measure RV using primarily 
the Mg-triplet (A ~ 5170 A) absorption feature; oth- 
ers use the infrared calcium triplet near 8500 A. Thus 
far we have used MMFS to measure 6415 RVs to a me- 
dian precision of ±2.0 km s~ 1 for 5180 stars in four MW 
dSphs. Of the measured stars, which span the magni- 
tude range 17 < V < 20.5, approximately 470 are mem- 
bers of the Carina dSph, 1900 are members of Fornax, 
990 are members of Sculptor, and 400 are members of 
Sextans. The remaining 1400 measured stars are likely 
dwarfs contributed by the MW foreground. In addition, 
the MMFS spectra are of sufficient quality that we can 
measure magnesium line strengths via a set of spectral 
indices that correlate with stellar parameters such as ef- 
fective temperature and surface gravity, and which trace 
chemical abundances of iron-peak and alpha elements. 
The spectral indices provide independent constraints on 



dSph membership and allow us to probe the chemical 
abundance distribution within a given dSph. 

Our focus in this paper is to provide a thorough de- 
scription of MMFS, the observations, data reduction 
procedure, derivation of measurement uncertainties, and 
comparison with previous velocity measurements. This 
material will be referenced in companion papers (Walker 
et al. in preparation) in which we present all data and 
analyses. 

2. THE MICHIGAN-MIKE FIBER SYSTEM 

In its standard slit observing mode the MIKE spectro- 
graph delivers high-resolution spectra spanning visible 
and near-infrared wavelengths. MIKE's first optical ele- 
ment is a coated glass dichroic that reflects/transmits 
incoming light into one of two independent channels. 
The "blue" channel has wavelength coverage 3200-5000 
A, while the "red" channel covers 4900-10000 A. The 
dual-beam design allows for independent optimization 
of throughput and dis persion charact e ristics over each 
range of spectrum (see Bernstein ct al. (2003) for a com- 
plete description of MIKE). 

For multi-object observations using Magellan+MMFS, 
MIKE is modified from its standard slit configuration 
to receive light from up to 256 fibers. In fiber mode 
MIKE is backed 2 m from the Nasmyth port, where it 
remains fixed with respect to gravity. Inserted at the port 
is a conical structure onto which an (interchangeable) 
aluminum plug plate is mounted. The plug plate holds 
fibers against a telecentrator lens that ensures the 1.4- 
arcsec fiber apertures accept light parallel to the optical 
axis and along the focal surface, regardless of position in 
the target field. A drum lens at the entrance to each fiber 
forms an image of the pupil on the the 175-^tm fiber core. 
The drum lens and the natural focal ratio degradation of 
the fiber convert the incoming beam from F/ll to F/3.5 
(more similar to MIKE's camera optics). Fibers run from 
the plug plate through a junction box, where they are 
bundled and sorted into one of two 128- fiber assemblies 
that enter the two channels of the spectrograph. Light 
losses are minimized by the short fiber length; from end 
to end, blue (red) fibers extend just 2.43 m (2.29 m). 

At the spectrograph end, the two 128-fiber assemblies 
replace MIKE's standard injection optics. Due to the 
spectrograph's fixed optical elements, blue and red chan- 
nels remain optimized for smaller and larger wavelengths, 
respectively. In order to fit spectra from all 128 fibers 
onto a single detector in each channel, narrow-band fil- 
ters are used to isolate a single order of spectrum for 
each target object. We use filters in both channels that 
isolate the spectral region 5130 — 5185 A. In late-type 
stellar spectra this region includes the prominent Mgl 
triplet (MgT) absorption feature and an assortment of 
Fel and Nil absorption lines. The MgT filter yields use- 
ful spectra through either MIKE channel, though the 
blue channel offers higher spectral resolution (~ 0.059 A 
pix -1 ; R ~ 25000) and greater throughput than does the 
red channel (~ 0.073 A pix" 1 ; R ~ 20000). MMFS thus 
has the ability to acquire high-resolution spectra from up 
to 256 objects simultaneously. The usable field of view 
subtends an angular diameter of 20 arcmin. 

3. OBSERVATIONS 
3.1. Target Selection 
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For our survey, MMFS target star candidates are se- 
lected from V,I imaging data taken prior to spectroscopic 
observing runs. We select MMFS target candidates from 
the dSph red giant branch (RGB) of the resulting color- 
magnitude diagram (CMD). We expect a fraction of the 
selected candidates to be foreground Milky Way dwarfs 
with magnitudes and colors placing them on or near 
the dSph RGB. The fraction of interlopers naturally in- 
creases for lines of sight along which the dSph surface 
density is relatively small. The efficiency with which 
our CMD-based selection identifies bona fide dSph red 
giants thus varies from galaxy to galaxy and with posi- 
tion within each galaxy, and foreground stars must sub- 
sequently be identified and removed from kinematic sam- 
ples. 

Details of the photometry used to select red giant can- 
didates in the Fornax and Sculptor dSphs are published 
elsewhere. Fornax candidates w ere selected using th e 
V,I photometric data described in I Walker et a l. (2006a). 
Sculpt or candidates were sele cted using the V,I photom- 
etry of IColeman et al.l (|2005l ). who generously provided 
their data set. 

3.1.1. Sextans 

For the Sextans dSph, red giant candidates were cho- 
sen based on V,I photometry taken with the MDM 2.4 
m Hiltner Telescope and 8k Mosaic CCD detector. Ob- 
servations took place the nights of 2004 February 11-17. 
Three of the seven nights were photometric; the remain- 
ing nights had mixed observing conditions. We observed 
25 Sextans fields, each with dimensions 24' x 24', form- 
ing a square array with adjacent field centers separated 
by 20'. We exposed for 360s in I and 540s in V. Stan- 
dard reduction steps included overscan, bias, and flat- 
field corrections using twilight sky fiats. We used the 
two-dimensional stella r photometry program DoPHOT 
(|Schechter et al.lll993l ) to identify stars and measure in- 
strumental magnitudes. We placed the photometry in 
the Kron-Cousins syste m (iBesselll 1 19761) using 50 pho- 
tometric standard stars (|LandoltJ ll9921 observed during 
the run. We used the 4' overlap between adjacent fields 
to remove photometric zero-point offsets. Formal error 
values returned by DoPHOT and multiple measurements 
from overlapping fields indicate the typical photometric 
accuracy is ±0.04 mag. 

3.1.2. Carina 

Over several MMFS observing runs, red giant candi- 
dates in the Carina dSph were chosen using three inde- 
pendent photometric data sets. Carina targets observed 
with MMFS in 2004 March were selected from B,R pho- 
tometry obtained with the CTIO 4m telescope with Mo- 
saic CCD detector during 1999 September. Carina tar- 
gets observed with MMFS in 2005 February were chosen 
in part from the same B,R mosaic photometry, but tar- 
gets in outer fields were selected from V,I photometry 
taken by Kaspar von Braun at the 2.4 m Swope Tele- 
scope (Las Campanas Observatory) in 2004 December. 
Carina targets observed with MMFS in 2006 March were 
selected from V,I photometry taken by Patrick Seitzer 
with the CTIO 0.9 m Schmidt Telescope in 2005 Decem- 
ber. We used DoPHOT to reduce the photometry in each 
case. 



In order to place the three photometric data sets for 
Carina on a common V,I system, first we used_l_22 photo- 
metric standard stars in the Carina field of Stetson (2000) 
to apply zero-point corrections to the Schmidt instru- 
mental V,I magnitudes. Residuals after applying these 
shifts have rms values 0.033 mag in V and 0.035 mag in I. 
We followed the same procedure in calibrating the 2004 
V,I photometry (residuals have rms 0.046 mag in V and 
0.048 mag in I). Next we used a piecewise linear transfor- 
mation to place the B,R mosaic photometry on the Stet- 
son V,I system. Coefficients were determined using 286 
Stetson standard stars detected in the mosaic photome- 
try. Over the color range 0.6 < V-I < 1.8, residuals have 
rms ~ 0.03 mag for the V and V-I transformations. We 
monitored for consistency by comparing the calibration 
results for stars detected in both the mosaic and Schmidt 
data. For 637 stars present in both data sets, flux devia- 
tions have mean values (V schm idt - V mosaic ) = 0.011 mag 
and (Ischmidt - Imosaic) = 0.003 mag, with rms values of 
0.072 mag and 0.088 mag, respectively. These values are 
reasonable considering the undersampling and blending 
in the Schmidt data. 

3.1.3. Globular Clusters 

For the purpose of calibrating spectral-line in- 
dex/abundance relations, we have also observed stellar 
targets in seven Galactic globular clusters. For clus- 
ters NGC 104 (47 Tucana), NGC 288, and NGC 7099 
(Messier 30) we chose targets on the basis of V,I pho- 
tometry obtained in 2005 July with the 40-inch telescope 
at Las Campanas Observatory, Chile. For clusters NGC 
5927, NGC 6121 (Messier 4), NGC 6171 (Messier 107), 
and NGC 6397 we selected targets using V,I photometry 
obtained in 2005 June with the Magellan/Baade 6.5m 
Telescope at Las Campanas. For all clusters, instrumen- 
tal magnitudes were measured with DoPHOT and then 
placed on the V,I system using stars in common with 
Stetson standard fields. For NGC 6171, which is absent 
from the Stetson catalog, we applied the same transfor- 
mations determined for NGC 5927, a cluster observed on 
the same night and at similar airmass. Direct compar- 
ison to Stetson stars indicates that the globular cluster 
photometry is accurate to within ~ 0.04 — 0.08 magni- 
tudes. 

3.2. CMDs and maps for dSph targets 

Figures[T]plots CMDs for Carina, Fornax, Sculptor and 
Sextans. Polygons enclose stars we consider to be red- 
giant candidates, and therefore eligible MMFS targets. 
For Carina, Sculptor and Sextans the selection region 
extends from near the tip of the red giant branch to the 
horizontal branch at V ~ 20.5. For Fornax, the greater 
stellar density allows us to increase efficiency by select- 
ing from only those red-giant candidates having V < 20. 
Figure [2] maps the locations of red giant candidates in 
each dSph. Circles in Figure [2] indicate MMFS fields 
that have been observed as of 2006 November. Initial 
MMFS runs gave priority to target-rich central fields, 
with goals of building sample size and evaluating instru- 
ment performance. Later runs gave preference to fields 
along the morphological major axis and/or at large ra- 
dius. The partial overlap of many adjacent fields allows 
us to quantify velocity variabilty using repeat measure- 
ments. 
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The 1.4" fiber apertures require that target coordinates 
be accurate to within ~ 0.2". DoPHOT's centroid algo- 
rithm returns the (x,y) CCD position of each identified 
star. We convert these to equatorial coordinates using 
the IRAF routines TFINDER and CCTRANS. We then 
tie the astrometry to the USNO-B1 system by identifying 
up to several hundred USNO stars per CCD frame. From 
the derived coordinates of the same stars in overlapping 
fields we find a typical rms scatter of 0.24". 

3.3. Spectroscopic Observing Procedure 

In a typical night we take 30-50 zero-second exposures 
just prior to evening twilight. During twilight we take 3- 
5 exposures of the scattered solar spectrum, to allow for 
wavelength calibration and provide a velocity zero point, 
followed by 3 x 60s quartz lamp exposures and a 1 x 300s 
Th-Ar arc-lamp exposure. 

Under optimal observing conditions we take 3 x 2400s 
exposures of target fields, immediately followed by 3 x 60s 
flat field (quartz) exposures and 1 x 300s Th-Ar. Other 
observations taken during the night include single-fiber, 
1 x 30s exposures of radial velocity standard stars, imme- 
diately followed by quartz and Th-Ar calibration expo- 
sures taken in the usual manner. The procedure allows 
for observation of 2-4 dSph fields per night. 

Table □ logs all MMFS observations taken for dSph 
and globular cluster fields as of 2006 November. The 
first two columns identify the field by galaxy name and 
field number (observed fields are mapped by field number 
in Figure [2|). Columns 3-7 in Table Q] list, respectively, 
the heliocentric Julian date at the midpoint of the first 
exposure, UT date at the midpoint of the first exposure, 
total exposure time, number of red giant candidates to 
which we assigned fibers, and the number of these for 
which the radial velocity was measured succesfully. The 
final column identifies repeat observations of fields for 
which we judged during the course of the run that the 
initial observation was insufficient. 

4. REDUCTION OF MMFS SPECTRA 

4.1. Initial Processing 

We extract MMFS spectra and measure radial veloc- 
ities using standard packages and subroutines available 
in the IRAF 4 astronomical software package. After 3x3 
binning of the 2k x4k detectors at readout, raw data 
frames from either detector are 683 x 1365 pix 2 , and 
the circular resolution element has FWHM ~ 3.5 binned 
pixels in the spectral direction. We reduce blue and 
red spectra independently and following identical pro- 
cedures. 

After applying overscan and bias corrections, we use 
the IMCOMBINE task to average the series of target 
frames for each field. Individual frames are scaled and 
weighted to account for differences in exposure time and 
airmass. Bad pixels, including most of those affected 
by cosmic rays, are rejected using a 4er clip around the 
median of the scaled frames. We then use the AP- 
SCATTER task to subtract scattered light from tar- 
get frames. We identify spectral apertures in the well- 
exposed quartz frame (averaged in the same manner as 

4 IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which is operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 



the target frames) associated with each target field, and 
then use the APALL task to extract one-dimensional 
spectra assuming the target, quartz, and Th-Ar frames 
associated with a given field have identical spectral aper- 
tures. 

4.2. Wavelength Calibration 

We determine the wavelength/pixel relation using ab- 
soprtion line features in the solar spectrum observed at 
evening twilight. In the observed region the solar spec- 
trum contains twice as many useful lines as the Th-Ar 
spectrum. The MMFS wavelength/pixel relation is suf- 
ficiently stable over an observing run that the benefit of 
using more lines outweighs disadvantages that arise be- 
cause twilight exposures are available only at dusk and 
dawn. To quantify any drift in the wavelength/pixel re- 
lation during the night we measure the relative redshifts 
of the twilight-calibrated Th-Ar spectra taken with each 
target field, and we compensate for these redshifts by 
applying low-order corrections to individual wavelength 
solutions. Because this method is a departure from the 
conventional Th-Ar calibration, we describe it in detail. 

In our twilight spectra we identify sol ar absorption 
lines using the NO AO FTS solar atlas (|Kurucz et al.l 
1984). Figure [3] displays an example twilight-solar spec- 
trum and labels lines used for wavelength calibration. We 
exclude several prominent lines (including the Mgl line at 
5167.327 A) because they are blended with nearby lines. 
The wavelength solution is determined from a 4* ft -order 
polynomial fit; residuals have mean rms ~ O.OlA (0.6 km 
s -1 ) for blue spectra and ~ 0.02A (1.2 km s -1 ) for red. 
For reference, typical residuals from fits to 9 emission 
lines in our Th-Ar spectra are larger by a factor of two. 

Next we apply initial wavelength solutions to all frames 
using the DISPCOR task. We assign all target, quartz, 
and Th-Ar frames the (aperture-dependent) wavelength 
solutions determined from the (temporally) nearest twi- 
light spectrum. We then use FXCOR to measure zero- 
point drift in the wavelength/pixel relation (during the 
night and/or over the course of the run) from shifts in the 
wavelengths of emission lines in the Th-Ar spectra. The 
measured zero-point shifts for Th-Ar frames taken during 
the observing run in 2005 February are shown as a func- 
tion of aperture number in Figure [4j We find that drift 
in the wavelength solution results in small (< 2 km s _1 ) 
zero-point offsets that depend linearly on aperture num- 
ber. We note that the zero-point offsets, and hence the 
wavelength/pixel relation, remain stable over the span of 
several nights. We correct for zero-point drift using the 
DOPCOR task to shift the wavelength solution in each 
aperture of the associated target and quartz frames by a 
value determined from the best linear fit to the zero-point 
offset vs. aperture data. 

Figure [5] displays solar radial velocity results measured 
from the several twilight exposures obtained during each 
MMFS observing run. Panels on the right side give ve- 
locity results for twilight-calibrated twilight frames. For 
comparison, panels on the left side give velocities mea- 
sured using conventional Th-Ar-calibration. We note the 
consistent velocity precision obtained with twilight cal- 
ibration. We suspect the degradation in the precision 
from Th-Ar calibration is related to the fact that the 
Th-Ar lamp system used exclusively by MMFS has op- 
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TABLE 1 

MMFS OBSERVATIONS OF DSPH AND GLOBULAR CLUSTER FIELDS AS OF 2006 NOVEMBER 
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Fig. 1. — V,I color-magnitude diagrams for the four dSphs observed with MMFS. Stars enclosed by polygons are considered red giant 
candidates and eligible MMFS targets 
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densely populated Fornax fields 4 and 9. 



erated at different voltages and with different power sup- 
plies over the six observing runs. The sample standard 
deviation of velocities measured using twilight calibra- 
tion has remained steady at ~ 0.8 km s _1 (blue) and 
~ 1 .6 km s _1 (red). Figure \5\ also provides evidence 
of small channel- and in some cases aperture-dependent 
zero-point offsets. The red channel tends to yield slightly 
more positive velocities than the blue channel. However, 
it is straightforward to measure and correct for these off- 
sets (see Section B~T]) . which are smaller than the typical 
velocity measurement error. 

4.3. Throughput Corrections and Sky Subtraction 

Because fibers differ in throughput by up to a factor of 
two, we create a response frame by dividing each aperture 
of the quartz frame by the average (over all apertures) 
quartz continuum. Dividing the target frame by the re- 
sponse frame corrects for variations in fiber throughput 



and pixel sensitivity. 

At this point we use the CONTINUUM task to remove 
from the target frames any remaining effects due to cos- 
mic rays and/or bad pixels and columns. We fit a 10 th - 
order legendre polynomial to the corrected spectrum in 
each aperture and replace with the function value any 
pixel values that deviate from the fit by more than 6a. 

Sky noise near the Mg-triplet is contributed mainly 
by scattered solar light. Left uncorrected, sky contam- 
ination produces a secondary velocity signal at km 
s . During observation of every dSph field we assign 
~ 32 science fibers (~ 16 per MIKE channel) to re- 
gions of blank sky. We use the SCOMBINE task to 
average the apertures corresponding to these sky spec- 
tra in the target frames. Prior to averaging we discard 
the three highest and two lowest values among the sky 
spectra at each pixel. Removal of the high pixel val- 
ues reduces the contribution from resolved background 




Fig. 3. — Example twilight solar spectrum used in wavelength 
calibration. Absorption lines used in the wavelength solution are 
identified. 



galaxies (|Wvse fc Gilmorel ll992). We subtract the aver- 
age sky spectrum from each aperture in the target frame. 

Black, solid lines in the left panels of Figures [6] and [7] 
display sky-subtracted spectra representing the range of 
S/N achieved during typical observations of dSph fields. 
Red, dotted lines in the left panels of either figure indi- 
cate the spectra prior to sky subtraction. The spectra in 
Figure [5] were obtained at new moon, while the spectra 
in Figure [7] were obtained at first-quarter moon. Prior 
to sky subtraction, spectra from the latter observation 
clearly are contaminated by solar absorption lines. Af- 
ter sky subtraction, residual solar lines typically deviate 
from the continuum level by an amount similar to the 
intrinsic continuum noise level. 

4.4. Final Processing 

In most MMFS observing runs there have been one or 
more fields for which we deemed the original observa- 
tion useful but insufficient (see Table 1). Reasons range 
from poor observing conditions to inadequate exposure 
time available prior to morning twilight. In such cases we 
re-observed the field using fiber assignments identical to 
those of the original observation. During the reduction 
steps already described, repeat observations — often sep- 
arated by several days from the original observation — are 
treated independently to allow for drift in the wavelength 
solution and different sky levels. At this point we com- 
bine all observations of the same field taken during the 
same run. First we use the DOPCOR task to correct 
the wavelength solutions of repeat observations for the 
changing component of Earth's motion (typical shifts are 
less than 1 km s _1 ) toward the target field. The SCOM- 
BINE task then averages spectra in corresponding aper- 
tures of the original and repeat target frames. Spectra 
in individual apertures are weighted by the mean pixel 
value prior to computing the aperture average. 

In the final processing step prior to velocity measure- 



ment, we remove the continuum shape of the spectra. 
The CONTINUUM task fits a 10 t,l -order legendre poly- 
nomial to the spectrum in each aperture of the target 
frames, then subtracts the function values from pixel val- 
ues to produce the final spectra. 

4.5. Signal-to- Noise Ratio 

The variance in our spectra has independent compo- 
nents contributed by source, sky background, and read 
noise. For the i th pixel, let N$ ,i represent the number 
of detected photons from the target star, let N B ,i repre- 
sent the number of detected photons from the sky back- 
ground, and let a\ be the noise associated with detec- 
tor readout (typically between 2-4 electrons per pixel for 
MIKE). In practice we estimate Ns,i from s individual 
sky spectra. Assuming the estimate of N B .i contributes 
a variance component N B ^/s, the signal-to-noise ratio 
(S/N) at the i th pixel is given by 



No 



V^ + N B ^ + N B ^/s + a 2 R (1) 

We estimate source and background photon fluxes from 
smooth fits to the the continua of our sky-subtracted and 
averaged sky spectra, respectively. At the i th pixel the 
continuum fit to the sky-subtracted spectrum has digi- 
tal unit (DU) value 5 n$,i — Ns,i/{RiG), where Ri is the 
value of the response frame (Section l4.3|) and G is the de- 
tector gain (between 0.5 and 1.5 electrons/DU over the 
several MMFS observing runs). The continuum fit to the 
averaged sky spectrum has DU value n Bj i = Nb,i/ (RiG). 
From the estimates Ns,i and Nb,% we calculate S/N at 
each pixel using Equation [TJ Figure [8] plots, for all spec- 
tra that eventually yield an acceptable velocity measure- 
ment, the distribution of the mean S/N calculated from 
pixels spanning the range 5141— 5177 A. The median S/N 
per pixel among dSph spectra acquired with MIKE's blue 
(red) channel is 3.7 (2.0) per pixel. Multiplying these val- 
ues by gives the mean S/N per resolution element. 

4.6. Measurement of Radial Velocities 

We measure radial velocity using the FXCOR task, 
which cross-correlates every continuum-subtracted spec- 
trum against a stellar template spectrum of known red- 
shift. For the template we use the sum of 428 spec- 
tra from bright, late-type stars, including 155 spectra of 
radial velocity standard stars and 273 spectra of stars 
from 7 Milky Way globular clusters. All spectra con- 
tributing to the template were obtained with echelle 
spectro graphs, either durin g a previous radial velocity 
study (jWalker et alj l2006al) or with MMFS as part of 
the present study, and all have spectral resolution equal 
to or higher than those of the target spectra. We co-add 
these individual spectra after shifting each to a common 
redshift. 

For cross-correlation, all spectra are re-binned linearly 
with log(A/A) at a common dispersion of 3.11 x 10~ 6 
pix -1 . Fourier transforms of the target and template 
spectra are filtered such that low-frequency continuum 
residuals and undersampled features at high frequencies 
do not contribute to the resulting cross-correlation func- 
tion (CCF). The right panels of Figures [6] and [7] display 

We follow the convention whereby q denotes the estimate of q. 
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Fig. 4. — Stability of the MMFS wavelength/pixel relation. Plotted for each multispectrum Th-Ar frame in the 2005 February MMFS 
run is the velocity shift measured via cross-correlation against the Th-Ar frame associated with a single twilight frame. Prior to cross- 
correlation, wavelengths in all Th-Ar frames were given identical wavelength/pixel solutions based on solar absorption lines identified in 
the twilight spectra. Blue/red points indicate blue/red MIKE detectors. Each panel indicates the heliocentric Julian date (—2.45 X 10 6 
days) at the midpoint of the Th-Ar exposure. 



CCFs calculated from the spectra shown in the left pan- 
els. The relative velocity between object and template is 
the velocity value at the center of a Gaussian profile fit to 
the tallest CCF peak (within a window of width 700 km 
s _1 , centered on the velocity of the target dSph). From 
the CCF, the known redshift of the template, the time 
of observation, and coordinates of the object, FXCOR 
returns the radial velocity in the solar rest frame. 

4.7. Velocity Zero Point 

The solar velocity measurements in Figure [5] show ev- 
idence of small residual offsets that depend on channel 
and aperture. We correct for these offsets by demanding 
that the many repeat velocity measurements obtained 
from solar twilight spectra yield mean velocities of km 
s _1 in each aperture in each run. We obtained sets of 
2 — 7 independent velocity measurements per aperture, 



per observing run, from solar twilight spectra (for a total 
of 6816 twilight spectra). We subtract the mean velocity 
of each such set from all velocity measurements (includ- 
ing those from twilight exposures, science targets and 
radial velocity standard stars) obtained using the same 
aperture in the same observing run; this places veloci- 
ties measured in a given run on a common zero point. 
The mean value subtracted from blue (red) apertures is 
0.48 km s _1 (1.07 km s _1 ), and the corrections have a 
standard deviation of 0.48 km s" 1 (0.78 km s" 1 ). The 
standard deviation from zero among corrected twilight 
RVs is 0.26 km s _1 (0.60 km s -1 ) for blue (red) aper- 
tures; we adopt this value as the baseline velocity error 
in the error analysis (Section [5j . 

Because the template spectrum is composed of late- 
type stellar spectra, velocities measured from solar spec- 
tra are susceptible to object /template mismatch. While 
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Fig. 5. — Comparison of the Th-Ar wavlength calibration 
method (left panels) with the twilight calibration method (right 
panels). For each of the MMFS observing runs, blue/red markers 
indicate (solar rest frame) radial velocities measured from twilight 
spectra in each aperture in blue/red channels. In 2005 October, 
Th-Ar exposures were not taken immediately after twilight expo- 
sures, precluding Th-Ar calibration of those twilight spectra. Due 
to the previous comparisons, we did not perform Th-Ar calibration 
of the 2006 October spectra. 

we consider the the zero-point corrections described 
above to be sufficient for removing channel-and aperture- 
dependent offsets within a given run, we next use dSph 
targets observed with the same channel in two or more 
runs to correct for variation in the velocity zero point be- 
tween observing runs. For such stars measured on blue 
(red) channels, Figure 191 (Figure [TU)) plots the difference 
between measured velocities as a function of the earlier 
measured velocity. Histograms indicate the distribution 
of measurement differences. In order to bring the cen- 
ter of each distribution to zero we apply constant offsets 
to all velocities measured in a given channel during a 
given run. Arbitrarily we choose the run/channel hav- 
ing the smallest mean velocities to serve as the absolute 
zero point. Table [2] summarizes the applied corrections. 
For each MMFS run listed in column 1, columns 2 and 
3 list the number of stars on blue and red channels that 
were observed in the same channel during at least one 
other run. Columns 3 and 4 list the blue and red offsets 
added to every velocity from the run. The maximum cor- 
rection is 1.7 km s _1 , slightly smaller than the median 
measurement error (Section [5]) . 

As a check on the absolute RV zero point we con- 
sider 58 velocity measurements of 22 RV standard 
stars (after applying the zero point shifts described 
above). Table ?? lists the individual MMFS velocity 
measurements for each standard, and for comparison 
lists previously publis h ed velocity measur e ments, from 
lOlszewski et al.1 (Il991h: lUdrv et all (jl999l ): IBeers et al.l 
(|2000h : iWalker etal (2006a|). The mean deviation be- 
tween the measured RV and the published value is 0.07 
km s _1 , which indicates the zero point of the MMFS ve- 
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locities is c onsistent wi t h RV standards. In the case of 
HD 45282, IBeers et al.l IpOOOh estimate a measurement 
error of ±10 km s _1 ; thus despite a deviation of 6 km s _1 , 
the MMFS measurement is consistent with the published 
value. Repeat MMFS measurements of RV standards 
show excellent agreement; the largest deviation among 
repeat measurements is 3.4 km s _1 , and 93% of repeat 
measurements differ by less than 3 km s^ 1 . 

4.8. Quality Control 

The right panels of Figures [H] and [7] demonstrate the 
range in quality of the CCF peaks we obtain from the 
MMFS spectra. In order to help eliminate poorly mea- 
sured velocities from the sample, we inspect all CCFs 
by eye and assign to each a pass/fail grade according to 
the following criteria. First, the tallest peak (within the 
search window centered on the dSph systemic velocity) 
in a satisfactory CCF must be unambiguous. CCFs in 
which the height of one or more secondary peaks reaches 
0.8 times that of the tallest peak are flagged as unsatis- 
factory. Second, the tallest peak in a satisfactory CCF 
must be symmetric about the measured velocity. Gaus- 
sian fits to asymmetric peaks can result in substantial 
velocity errors; we therefore flag as unsatisfactory those 
CCFs having significant residuals with respect to the 
Gaussian fit. According to these criteria, the CCFs in 
the top two panels of both Figure [6] and Figure [7] are 
considered unsatisfactory. 

A second, quantitative quality control filter is pro- 
vided by the Tonry-Davis value, Rtd (calculated by 
FXCOR), which indicates the height of the tallest CCF 
peak relative to the average CCF peak (|Tonrv fc David 
fl979h . Fi gure [11] plots, separately for all the dSph stellar 
and blank-sky spectra measured in blue and red chan- 
nels, the distributions of Rtd associated with CCFs 
judged by eye to be satisfactory (solid histograms) and 
unsatisfactory (dotted histograms). We expect blank- 
sky spectra to yield poor CCFs; this is generally the 
case, as indicated by the small number of sky spectra 
judged by eye to give satisfactory CCFs and the rela- 
tively low Rtd values associated with these measure- 
ments. Among stellar spectra, qualitatively satisfactory 
CCFs have 2.8 < R T d < 40, with median Rtd = 9.12. 
We define a critical value, Rtd,c, to be the value of Rtd 
at which, absent other information, the probability that 
a stellar CCF is satisfactory equals the probability that it 
is unsatisfactory. For spectra obtained with MIKE's blue 
channel this occurs at Rtd,c = 4.0; for the red channel 
Rtd,c = 4.5. On the basis of Figure [TT1 we allow a given 
stellar velocity measurement into the final sample only if 
it is derived from a CCF that 1) is judged by eye to be 
satisfactory, and 2) has Rtd > Rtd.c- 

5. VELOCITY MEASUREMENT ERRORS 
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Fig. 6. — Processed spectra and corresponding radial velocity cross-correlation functions from observations during a dark night. Solid 
black lines in the left-hand panels display examples of fully processed spectra from a Sculptor target field. Red lines represent each 
spectrum prior to sky subtraction. Right-hand panels show each spectrum's cross-correlation against our radial velocity template and list 
the Tonry-Davis R value measured from the CCF. 



5.1. Gaussian Errors 

For the majority of stars in our sample we have ob- 
tained only a single velocity measurement and direct es- 
timation of the measurement error is impossible. Instead 
we use the many existing cases of independent repeat 
measurements to relate the errors to quantifiable features 
of individual spectra. Let Vij be the i th of nj indepen- 
dent radial velocity measurements obtained for the j th of 
N stars in the data set, and suppose the j th star has true 
radial velocity V*j . We adopt a model in which the veloc- 
ity measurement error, Vij — V*j, is the sum of two inde- 
pendent Gaussian components. First, noise in the stellar 
spectrum gives rise to noise in the cross-correlation func- 
tion. We assume the error component due to CCF noise 
follows a Gaussian distribution with variance Occf- The 
finite spectral resolution of the instrument contributes a 
second error component that is independent of signal- 
to-noise ratio and defines the maximum precision of the 



instrument. We assume this "baseline" error follows a 
Gaussian distribution with constant variance, Uq. The 
error in a given velocity measurement is then 



Vij - V*j = (o- C CF,ij + a o) 



(2) 



The values e^, which account for the inherent random- 
ness of measurement, follow the standard normal (Gaus- 
sian with mean zero and unit variance) distribution. The 
resulting error bars are symmetric about Vij, and given 
by ±oy,y, where oy,y = (occfm + ° o) 1/2 



We model accF as in Wal ker et alJ (|2006al) : 



{1 + RTD^jY 



(3) 



Except for the introdu ction of the par a meter x, this is 
essentially the model of iTonrv fc Davisl £1979), in which 
a C cF tx (1 + RtdY 1 - 
After substituting for cjccf in Equation [2j taking the 
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Fig. 7. — Panels plot the same relations as in Figure[6] but for example spectra from a Sextans field observed in the presence of significant 
moonlight. 



base-ten logarithm gives 



log[(Vy - V*i) 2 ] = log 



{I + Rtd, 



\2x 



log[e 



(4) 

For normally distributed , the values ef 3 follow the x 
distribution (one degree of freedom) . From Monte Carlo 
simulations, log[e? ] has mean value (log[e? ]) = —0.55. 



mizes the sum 

N ™j , 

EE 1 * 

j=i t=i v 



-v.. 



log 



a£ 1-0.55 



If we define cL. 



log[e? ] + 0.55, then 
log 



\og[(V l3 -V^ 2 } = log[ (1 + flTDt . )2g +^j +*«-0.55, 

(5) 

and (Sij) = 0. For the baseline variance we adopt the 
standard deviation of velocities measured from twilight 
spectra. The CCFs from these spectra generally have 
Rtd > 30, a regime in which accF is negligible. From 
Sectionl^Tl a = 0.26 km s" 1 (cr = 0.60 km s _1 ) for blue 
(red) channels. The parameter pair {x, a} can then be 
estimated, via nonlinear regression, as that which mini- 



( ) 

We estimate {cc, a} separately for blue and red chan- 
nels. We consider only dSph candidate stars with multi- 
ple, independent velocity measurements that pass quality 
control filters fSection l4.8[) . For the blue (red) channel we 
use 1249 (561) measurements of 583 (257) stars. We re- 
place the unknown V*j with an estimate computed from 

the weighted mean, V*j — Y^%=1 w ijVij I Ym=i w v '•■ For 
weights we use Wij — {ctpxij + °o) _1 j where <JFX,ij is 
the initial velocity error returned by FXCOR (to first 
order, <JFX,ij is proportional to <JccF,ij)- This substitu- 
tion introduces a bias, as it ignores the variance in the 
estimat e V*j . However, analysis of repeat measurements 
(Section l5.3|) demonstrates that resulting error estimates 
are reasonable. 
The top panels in Figure [12] plot (for blue and red 
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TABLE 3 

MMFS Velocities of Radial Velocity Standards 



ID 


HJD 


^published 


V 


^ ^published 




(-2.45 xlO 6 days) 


(km S" 1 ) 


(km s — 1 ) 


(km S — 1 ) 


HD 83212 


3088.6241 


110.0 ± l.O 3 


107.4 ±0.4 


-2.6 ± 1.1 




3088.6276 


110.0 ± 1.0 


107.5 ± 0.4 


-2.5 ± 1.1 




3088.6321 


linn 1 1 /~\ 

110.0 ± 1.0 


109.2 ± 0.7 


—0.8 ± 1.2 




3088.6363 


nn n _J_ 1 n 
110.0 ± 1.0 


1 nn q _l_ n *? 
109.3 ± 0. / 


—0.7 ± 1.2 


HD 103545 


3088.6441 


177.0 ± 1.5 4 


178.9 ± 3.3 


1.9 ± 3.6 




3088.6491 


177.0 ± 1.5 


176.9 ± 1.3 


-0.1 ± 2.0 




3090.6633 


177.0 ± 1.5 


177.6 ± 1.2 


0.6 ± 1.9 




oUyU.DDDD 


1 77 n _L 1 c; 

l / i.u ±1.0 


1 7fl c; _L Q A 

1 /y.o ± 0.4 


2.5 ± 3.7 




Q/LOQ R8.49 

o'±uy.oo i iz 


I 77 n _l 1 c 

II / .u m 1.0 


1 77 9 _|_ 1 c: 
1 / ( .Z X l.O 


2 ± 2 1 




3417 8769 

-LI . C? 1 UC/ 


177 ± 1 5 


176 9 ± 1 2 


—0 1 ± 19 




3418.8698 


177^0 ± L5 


177^8 ± l!4 


ois ± 2!o 


HD 45282 


3089 4Q00 


Qf)i + 1 3 


307 14-24 


6.1 ± 10.3 


xlJJ oUl I U 


oUoy.DOOo 


U.O ± U.z 


n ^ 1 n 7 
U.O ± U. / 


n n _l_ n 
U.U ± U.o 




oUoy.uoo4 


u.o ± u.z 


i a J- n q 
— l.O ± U.v5 


— Z.l ± U.4 




QfiOO fiQ 1 A 


u.o ± u.z 


— z.o ± U.o 


1 n a 
— 0.0 ± U.4 




oUyz.oooo 


0.5 ± 0.2 


— 0.1 i 0.7 


—0.6 ± 0.8 


HD 92588 


3089.6641 


42.5 ± 0.3 2 


40.8 ± 0.4 


— 1.8 ± 0.5 




ouoy . uuoo 


40 _l n 3 


42.1 i 0.7 


—0.4 ± 0.8 




3091.6315 


42.5 ± 0.3 


41^9 ± 0.7 


—0.Q ± 0.8 




3091 6343 


42 5 + 3 


41.1 ± 0.3 


— 1.4 ± 0.4 


TT"P) 1 OA Q c; Q 
xlJJ iZ4o0o 


QflOn QC\RA 

oUyU.oUo4 


ozo.u ± i.u 


QOQ P 1 O 7 


— 1.4 ± z.y 




ouyu.oii i 


OZO.U HI J. .u 


qoQ 4-07 
ozo.u zc u. / 


2 ± 1 2 


TTTl 1 ^7A^7 


QflQfl SI 1 ^ 

ouyu.oiio 


17S4-n 
K.O I U.o 


1 « 1 _i_ n ^ 

ID. 1 I U.O 


1 1 4- n A 




3090 81 34 


17.8 ± 0.3 




— 1.2 ± 0.8 


xliJ iloUOO 


oUyl.D4ol 


— 1U1.U ± l.U 


— lUl.o ± U.4 


— U.o ± 1.1 




ouyi.U'± / i 


1 ni 4-1 n 

— 1U1.U HI l.U 


1 09 9 4- n Q 
— iuz.z ^ u.y 


1 9-1-1 Q 
— ±.z in ±.0 






1 ni + 1 


— 103 2 + 04 


—2.2 ± 1.1 




3417 8831 


— 101 0-1-1 


— 102 7 + 04 


— 1^7 ± LI 




3801 8875 


— 101 + 1 


— 100 2 + 4 


8 ± 11 




3801.8891 


-101.0 ± 1.0 


-100.1 ± 0.4 


0^9 ±1/1 


HD 176047 

±1.1— s _1_ p l^Urt 1 


3287 5020 


—42.5 ± 0.2 2 


—41.8 ± 0.3 


0.7 ± 0.4 


TTTl 1 QAQS^ 

lLLj lyoyoo 


Q007 c:i 1x7 
0Z0 / .010 / 


— y . 1 x u.o 


s -t- n ^ 


U.o III U.4 




4018 5032 


—9 1 ± 3 


— 9 3 ± 3 


—0 2 ± 4 




4027 4955 


-9+ ± 0^3 


-9^8 ± 0^3 


-0^7 ± 0^4 




4027.5121 


-9.1 ± 0.3 


-7.0 ±0.7 


2.1 ±0.7 


HD 219509 


3287.5246 


67.5 ± 0.5 2 


66.3 ± 1.0 


-1.2 ± 1.1 


CPD-432527 


3287 8984 


19.7 ± 0.9 2 


20 4 + 3 


0.7 ± 0.9 




341 51 82 


19.7 ±09 


17.9 ± 0.3 


— L8 ± 0^9 




3665 8825 


19.7 ± 0^9 


19 3 + 3 

J-C.O _1_ V J , • J 


— 0^4 ± 0^9 




ozoo.oy4o 


4u.o IE u.z 


1 .0 IE U.O 


1.1 III u.o 




3289 8915 


40.5 ±02 


43.2 ± 0.3 


2.7 ± 0.3 




3289.8939 


40.5 ± 0.2 


43.4 ±0.3 


2.9 ±0.3 




3666.8791 


40.5 ± 0.2 


39.9 ±0.3 


-0.6 ±0.3 


SAO 201636 


3408.8782 


264.3 ± 1.4 4 


264.8 ± 1.0 


0.5 ± 1.7 




3409.8670 


264.3 ± 1.4 


264.2 ± 0.7 


-0.1 ± 1.6 


HD 83516 


3409.8590 


43.5 ±0.2 2 


41.8 ±0.3 


-1.7±0.3 




3412.9035 


43.5 ± 0.2 


42.0 ±0.3 


-1.5 ±0.3 


HD 111417 


3410.8967 


-19.1 ±0.2 2 


-19.1 ±0.6 


0.0 ±0.7 




3803.8877 


-19.1 ± 0.2 


-18.7 ±0.6 


0.4 ±0.7 


HD 43880 


3411.5240 


43.6 ± 2.4 1 


40.6 ±0.3 


-2.9 ±2.4 




3412.5257 


43.6 ± 2.4 


40.3 ±0.3 


-3.3 ±2.4 




3801.5077 


43.6 ± 2.4 


43.3 ±0.3 


-0.3 ±2.4 




3801.5090 


43.6 ± 2.4 


43.2 ±0.3 


-0.4 ±2.4 


HD 93529 


3411.9006 


143.0 ± 1.0 3 


144.3 ±0.3 


1.3 ± 1.1 


HD 23214 


3415.5220 


-5.1 ±0.9 4 


-6.8 ±0.3 


-1.7 ± 1.0 


HD 21581 


3416.5260 


151.3 ± 1.3 4 


151.7 ±0.5 


0.4 ± 1.4 


SAO 217998 


3417.5155 


19.0 ± 1.6 4 


17.9 ±0.3 


-1.1 ± 1.6 


HD 223311 


4018.5258 


-20.2 ±0.3 2 


-20.0 ± 0.4 


0.2 ±0.5 



1 IQlszewski et al.l II1991I1 2 lUdrv et al.l IH999D 3 iBeers et alj (|2000f) 4 Imiker et al.l ll2006af) 



channels) the nonlinear regression data in the log[(Vy — 

V*j) 2 ] versus (1 + RTD,ij) plane. Overplotted is the best- 
fit curve. We note the presence of outliers that have 
\V~ij — V*j \ up to ~ 250 km s _1 (middle panels of Figure 
[12p . The broad velocity distribution of these outliers 
suggests an error component not considered in the model 
specified by Equation^ We address non-Gaussian error 
in Section [5\2l 

Table |4] summarizes the estimated parameters of the 
Gaussian error model. For the blue (red) channel, the 
best-fit curve to the regression data corresponds to the 
parameter pair {x, a} — {1.69, 81.0 km s -1 } ({1.65, 107.9 
km s -1 }). Using these values, we calculate error bars 
±a y^j as a function of Rtd (bottom panels of Fig- 
ure [T2]) . Velocities measured from minimally acceptable 
CCFs (those with Rtd — Rtd,c), have <Jv,ij = 5.3 km 
s _1 (blue) and oy,ij = 6.5 km s _1 (red). Among all 



acceptable measurements, median errors are ay.ij = 1-6 
km s^ 1 (blue) and uv.ij = 2.4 km s _1 (red). 

5.2. Non-Gaussian Errors 

The adopted error model assumes the measurement er- 
ror is the sum of Gaussian components. Figure [T2l shows 
that the several outliers from the main distribution of 
regression data (top panels) correspond to outliers in the 

Vij — V*j distribution (middle panels). The distribution 

of Vij — V*j among these outliers indicates an additional, 
non-Gaussian error component not considered in Equa- 
tion 2. For stars with repeat observations it is straight- 
forward to identify and discard affected measurements; 
this is not possible for stars with only a single measure- 
ment. However, we show that only a negligible number 
of measurements affected by large non-Gaussian errors 
are likely to be included in dSph samples. 
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TABLE 4 

Summary of Parameters for Gaussian Error Model 



Channel 


ttTD.c 


^ repeats 


CO 


X 


a 


^V, median 










(km s- 1 ) 




(km s _1 ) 


(km S _1 ) 


(km s _1 ) 


Blue 


4.0 


1249 


0.26 


1.69 


81.0 


1.6 


5.3 


Red 


4.5 


561 


0.60 


1.65 


107.9 


2.4 


6.5 





il !0 )) « 50 

S/N 




Fig. 8. — Distributions of the mean signal-to-noise ratio (S/N) 
per pixel, calculated over the wavelength range 5141 — 5177 A. 
Distributions are plotted separately for dSph, globular cluster, RV 
standard and solar twilight spectra. In each panel the thin solid 
(dotted) lines give distributions for spectra obtained with MIKE's 
blue (red) channels, while the thick solid line is the combined dis- 
tribution. The resolution element is 3.5 pixels and the dispersion 
is 0.08 A/pix (0.12 A/pix) for blue (red) channels. 
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Fig. 9. — Distribution of measurement deviations for stars mea- 
sured with the blue channel in multiple observing runs. Left-hand 
sub-panels plot velocity deviation as a function of the (chronolog- 
ically) first measurement; histograms in the right-hand subpanels 
give the distribution of AV. We observed a given star in up to 
three different runs. Filled circles and solid histograms result from 
comparison with the earlier of the two other runs; open circles and 
dotted histograms result from comparison with the latter of the 
two other runs. 



Potential sources of non-Gaussian errors include astro- 
physical velocity variability (e.g., binary stars), but the 
most plausible origin of large Vij — V*j outliers is the se- 
lection of a false peak from a noisy CCF. The outlying 
data points occur preferentially at small Rtd, a regime 
in which CCFs are relatively noisy and odds of selecting 
a false peak increase. Adjacent peaks in CCFs measured 
from our spectra are separated typically by 50 — 100 km 
s" 1 . Assuming false peaks occur at no preferred velocity 
within the CCF, velocities derived from false peaks will 
result in large deviations that follow a uniform velocity 
distribution. This is consistent with the observed dis- 
tribution of Vij — V*j among outliers (middle panels in 



Figure [T3 

To what extent does false peak selection affect our data 
set? Among the dSph candidate stars with multiple mea- 
surements, we consider cases of false peak selection to be 
associated with measurements having \Vij — V*j\ > 30 
km s _1 and Rtd < 6; we do not consider Vij — V*j 
outliers with Rtd > 6 to be false-peak measurements 
because all such cases can be attributed to a poor es- 
timate of V,j due to an associated repeat measurement 
having Rtd < 6. For the blue (red) channel, 22 of 1249 
(6 of 561) repeat measurements that pass both of the 
initial quality control filters (Section |4.8|) meet the false- 
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Fig. 10. — Same as Figure [9] but for stars observed with the red 
channel in multiple runs. 
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Fig. 11. — Distributions of Rtd associated with measured cross- 
correlation functions. Panels indicate whether the distribution cor- 
responds to stellar or blank-sky spectra obtained in MIKE's blue or 
red channel. In each panel, the solid/dotted histogram is the dis- 
tribution of Rtd associated with CCFs judged by eye inspection 
to be satisfactory /unsatisfactory. 



E 

i 






Fig. 12. — Top: Data from repeat velocity measurements, 
used for estimating error model parameter pair {x, a}, with best-fit 
curve overplotted. Left/right panels show data for blue/red chan- 
nels. Middle: Deviation of measured velocity from estimated true 
velocity. Bottom: Size of error bar resulting from parameter pair 
{x, a} that yields the best-fit curve. In all panels, dotted, vertical 
lines indicate the minimally acceptable value of Rtd (section 14. 81 . 



peak criterion. Combining results from both channels 
yields 486 additional repeat measurements that pass the 
quality control filters; ten of these meet the false peak 
criterion. Extrapolating to the entire sample, we expect 
fewer than 2% of measurements passing quality control 
filters to be derived from false CCF peaks. Cases that 
are identifiable via comparison to repeat measurements 
are removed from the sample, but we expect ~ 70 cases 
of false peak selection to remain among the 4175 stars for 
which we obtained only a single measurement. However, 
if the distribution of unidentified false-peak velocities is 
uniform over a range spanning at least 400 km s _1 , as 
Figure [T2l suggests, we expect only 10—15% to fall within 
the range of a typical dSph velocity distribution. This 
results in only a handful of false-peak velocities spread 
among the four dSph member samples. This level of con- 
tamination is below that expected to be introduced by 
cases of ambiguous dSph membership. 

5.3. Repeat Measurements 

After removing identified cases of false CCF peaks, 
the final MMFS data set contains repeat velocity mea- 
surements for 1011 dSph target stars. There are 2246 
independent measurements of these stars, including up 
to five measurements for some stars. If we calculate 



the weighted mean velocity for the j th 
-n) — - — — , using weights wu = o\ 



-2 
'V,ij> 



star as Vj — 
then the Val- 



ues Vij - 
variance 



Vj should follow a Gaussian distribution with 



4k. 3 = Var(Vij 



wu 1 K ' K 2 t—* 



where K = 



V) = — I J - — ) + — > w k j, 

(7) 

Figure [T3] plots the distribution of 
(Vij ~ Vj)l a i±v,ij for all dSph targets with repeat velocity 
measurements. A Gaussian fit to the empirical distribu- 
tion has standard deviation a — 1.28. The excess of this 
value with respect to the nominal value of one is due in 
large part to the presence of outliers and is not surprising 
given the likely presence of binaries in the sample. 

6. COMPARISON WITH PREVIOUS WORK 

The MMFS sample has stars in common with 
several prev iously pub lish ed RV data sets — 
Carina: IMateo et al.l (|1993f) ; iMaiewski et al 



1991ft: I Walker et all J2006aft 



(|2005ft : iMunoz et all (12006ft. F ornax: IMateo et a! 



Sc ulptor 

lArmandroff fc Da Costal (11986ft ; iQueloz et all (1 1995ft 



Sextans: iDa Costa et all (11991ft: iSuntzeff et all (jl993li : 
lHargreaves et al.l (|1994ft: iKlevna et all (12004ft. Among 
previo us st udies, only IMateo et alj (|1991ft . iQueloz et al.l 
(|1995ft . and lWalker et alj (|2006aft measure RV using the 
magnesium triplet. All others use the infrared calcium 
triplet (CaT) near 8500 A. The MMFS sample offers the 
first opportunity to compare directly large numbers of 
velocity measurements obtained using both techniques. 

To identify MMFS stars in common with other sam- 
ples, we searched published data sets for coordinate 
matches, tolerating offsets of up to 0.2" — 1.5" depending 
on the sample. For each detected matching pair, Figure 
[14] plots the velocity deviation AV = Vmmfs — V^ther as 
a function of Vmmfs- Figure [T5l plots the same quan- 
tities over a window of width 50 km s _1 , centered on 
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300 




(V-V)/<r lv 

Fig. 13. — For all stars with repeat MMFS measurements, the 
distribution oiV — V, normalized by the estimated error. Overplot- 
ted is the best-fitting Gaussian distribution (solid line), which is 
nearly identical to the (vertically scaled) normal distribution with 
unit variance (dotted line). 



the systemic velocity of the dSph. Also shown in both 
figures are internal comparisons using stars with repeat 
MMFS measurements. In each panel of Figures [14] and 
1151 the marker type indicates the relative quality of the 
RV measurement. Filled squares denote stars for which 
both the MMFS and comparison measurements have 
quoted uncertainties less than the median uncertainty 
in their respective sample. Open squares/triangles de- 
note stars for which the MMFS measurement has un- 
certainty less/greater than the median MMFS uncer- 
tainty and the comparison measurement has uncertainty 
greater/less than the median from its sample. Crosses 
denote stars for which both the MMFS and compari- 
son measurements have uncertainties greater than the re- 
spective sample median. We find that the best measure- 
ments (solid squares) indeed provide the closest agree- 
ment between samples. Figure [16] plots for each compar- 
ison sample the distribution of AV/ \J c^[ MFS + Co thcr - 

Table [5] summarizes the comparisons. For each com- 
parison sample, columns 3-6 list the primary absorption 
features used, the spectral resolution (if available), the 
number of stars in common with the MMFS sample, 
and the number of MMFS measurements for the com- 
mon stars. Columns 7-8 list the mean velocity offset and 
standard deviation calculated from all matching pairs. 
Because these values are susceptible to the effects of even 
a single outlier (which may be due to a mis-identified 
match, a binary star, or a false CCF peak) , columns 9- 
10 list biweight estimates (jBeers et alJll990( ) of the mean 
offset and scatter. While the biweight is resistant to out- 
liers, it is known to perform poorly when the number of 
matches is small (N < 10). 

We find generally acce p table agreement wi t h the 
resul ts of iDa Costa et al.l til991l); iMateo et al l (11991L 
119931) : ISuntzeff et alj (|1993f k~ lHargreaves et all (|1994D ; 



iQueloz et all (|1995t ); iMaiewski etTal! (l2005fh Despite 
large a ppa rent deviations, we incl ude the IDa Costa et all 
(|1991h and lMaiewski et all (|2005l) samples in this list be- 
cause both studies report large measurement uncertain- 
ties (10-20 km s _1 and 4-10 km s _1 , respectively). The 
somewhat l arge offset and stan dard deviation with re- 
spect to the lQueloz et al.l (|1995l) Sculptor sample are due 
to a single star; the remaining five matching stars show 
excellent velocity agreement. The single MMFS star in 
common with the lArmandroff fc Da Costal (|1986f) Sculp- 
tor sample happens to be the same star giving the most 
deviant result in the Queloz et al. comparison; notice 
that for this star the Armandroff et al. and Queloz et al. 
measurements agree. Unless the MMFS measurement is 
in error, it is likely that this star is a binary, consistent 
with the conclusion of Queloz et al. that Sculptor has 
a binary fraction of ~ 20%. Other than this, we can 
conclude little from the s ingle star in common with the 
lArmandroff fc Da Costal (|1986f ) sample. Similarly , the 
three stars in common with the Mateo et al.l (1991) For- 
nax sample offer little basis for comparison. 

The comparison to the Fornax sample of IWalker et al.l 
(2006aj) shows evidence for a zero-point offset of ~ 2 
km s , and the scatter is somewhat larger than ex- 
pected from quoted measurement errors. Here we can 
rule out mis-identified coordinate matches, as MMFS 
targets are c h osen f rom the same target list used by 
IWalker et al.l (|2006aD . While it remains possible that 
some disagreement is due to binary stars, we note that 
the most deviant measurement s correspond to stars for 
which the Wal ker et alj (|2006af ) measurement has higher 
uncertainty than the median for that sample . If the mea- 
surement errors for the lWalker et al.l (|2006af) stars are 1.3 
times larger than quoted, the distribution of AV becomes 
that which we expect for Gaussian errors. We note that 
the erro r model given by Equ ation [^improves upon that 
used in IWalker et al.l (|2006aD . which does not take the 
baseli ne variance into acco unt. 

The lKlevna et alj (|2004l ) and lMunoz et alj (|2006l ) sam- 
ples offer relatively large numbers of common stars (45 
and 229, respectively) with which to compare MMFS ve- 
locities. The Munoz et al. sample consists of two data 
sets acquired independently using Magellan+MIKE (in 
slit mode; 11 MMFS stars) and VLT+GIRAFFE (fibers; 
218 MMFS stars); we consider each sample separately. 
The most striking feature we find in these comparisons 
is the apparent correlation of AV^ with Vmmfs- With 
respect to the VLT+GIRAFFE data, this manifests as 
a tendency toward (Vmmfs — Wlt) < km s _1 when 
I^mmfs < 75 km s _1 (Figure [14|). We note that the stars 
with largest VLT uncertainties and smallest MMFS un- 
certainties (open squares) are most responsible for this 
trend. Moreover, for the several stars that have incom- 
patible MMFS and VLT velocities and have multiple 
MMFS measurements, we find that the MMFS measure- 
ments show good internal agreement. 

In compari son to the Mu n oz et al. MIKE data for Ca- 
rina and the iKlevna et alJ (|2004T ) data for Sextans, we 
find apparent correlations between AV and Vmmfs over 
the velocity distribution characteristic of the dSph mem- 
ber stars. Dotted lines in each panel of Figure [15] give 
least-squares linear fits to the plotted points. Columns 
11-12 in Table[5]give the slope of each line and the proba- 
bility that such a slope would be observed in the absence 
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TABLE 5 

Comparison of RV data with previous results 



Galaxy 


Author 


A 


.Resolution 




N v 


1 A 1 I \ 

(A.V) 
(km s ) 


ft — 1 \ 

(km s ) 


AV BW 
/i — i \ 
(km s 1 j 


CT BW 
(km s 1 ) 


slope 


P 


— 7"i ' 

Carina 


Mateo ct al. (1993) 


— in — ttt — 

MgT 


21000 


5 


6 


-2.34 


3.54 


-0.45 


3.75 


0.27 


0.31200 


Carina 


Maiewski et al772005l 


CaT 


2600,7600 a 


17 


20 


40.54 


123.32 


-2.04 


12.51 


-0.23 


0.39328 


Carina 


Munoz et al. (20061 a 
Munoz et al. g006) c 


CaT 


19000 


11 


17 


2.28 


7.94 


0.62 


5.06 


0.49 


0.01460 


Carina 


CaT 


6500 


218 


252 


-0.98 


26.88 


0.51 


6.05 


0.13 


0.15372 


Sextans 


Da Costa et al~( 1991) 


CaT 


2600 


6 


14 


-5.58 


14.88 


-2.41 


20.71 


1.15 


0.08198 


Sextans 


Suntzefl et al. ( 1993 ) 


CaT 


2800 


4 


14 


-5.53 


6.35 


-1.86 


7.50 


0.018 


0.88844 


Sextans 


Hargrcavcs ct al. (1994) 


CaT 




11 


32 


-0.42 


5.26 


-0.61 


4.64 


0.050 


0.68364 


Sextans 


Klevna et al. (2004) 


CaT 




15 


92 


-6.84 


24.80 


-0.80 


5.46 


0.34 


0.00178 


Fornax 


Mateo et al. (1991) 


MgT 


21000 


3 


3 


2.41 


2.82 


0.40 


1.40 


0.37 


0.14845 


Fornax 


Walker et al. (2006a) 


MgT 


21000 


58 


72 


3.33 


7.76 


1.62 


5.23 


0.094 


0.07913 


Sculptor 


^^mandrolii^^^^osta (1986) 


CaT 




1 


1 


13.32 












Sculptor 


Queloz et al. (1995) 


MgT 


16000 


5 


6 


2.94 


6.56 


0.42 


2.62 


0.90 


0.02979 



Instrumental setup differed over two observing runs; sec Maiewski ct al. (2005) for details. ^ using Magcllan+MIKE C using VLT+GIRAFFE 



of a correlation (determined via 100000 Monte Carlo re- 
alizations that assume the same Vmmfs values and draw 
AV randomly from a Gaussian distribution with vari- 
ance equal to the standard deviation of AV for points 
plotted in Figure I15[) . The apparent trends observed 
with respect to the Munoz et al. (MIKE) and Kleyna 
et al. comparisons have slopes d( AV)/dVMMFS = 0.49 
and 0.34, respectively, and these are significant at the 
98.5% 99.8% levels. Over the range of velocities con- 
sistent with dSph membership, these slopes amount to 
systematic deviations of up to 10 km s _1 . 

That we observe such behavior in the largest available 
comparison samples suggests there may be systematic 
differences between survey data. At this point we can 
only speculate regarding the source of such a system- 
atic trend and whether it might originate in the MMFS 
or the two previously published samples. Perhaps the 
most significant observational difference between ours 
and the other large dSph surveys is that we observe at 
the MgT rather than the CaT. Whereas late- type stars 
are brighter in the infrared, the primary benefit we enjoy 
by using the MgT is the relatively clean sky. In contrast, 
the CaT region is prone to contamination by atm ospheric 
-OH emission lines (see lOsterbrock et al.lll996h . At red- 
shifts characteristic of Carina and Sextans (~ 225 km 
s" 1 ), strong s ky lines fall di r ectly on top of lines 1 and 
2 of the CaT. IKlevna et all (|2004D note that, when us- 
ing CaT lines individually to measure velocities, lines 1 
and 2 tend to give spurious Sextans-like velocities of 233 
and 223 km s , respectively, due to contribution from 
sky residuals to the cro ss-correlation f unctio n (see their 
Figure 1). In the end, IKlevna et all (|2004D use lines 2 
and 3 jointly for their velocity measurements. If the in- 
clusion of line 2 tends to "pull" their velocities toward 
223 km s _1 , this could account for the observed gradient 
in Figure [T51 Munoz et al. (w/ MIKE) measure Carina 
velocities using ~ 15 absorption features, including all 
three CaT lines. Because Carina's radial velocity in the 
solar rest frame is nearly identical to that of Sextans, one 
might expect a similar effect from sky residuals. Indeed 
we observe in comparison to the Munoz et al. (w/MIKE) 
sample a slope similar in magnitude and significance to 
those obtained in the comparison to the Kleyna et al. 
sample. 

7. SPECTRAL LINE INDICES 

Our MMFS observing program was designed chiefly to 
measure stellar velocities at low S/N. While our spectra 



are not suitable for measuring equivalent widths, nearly 
all are of sufficient quality for estimating line strengths 
via spectral indices. Let S(X) represent the observed 
spectrum and suppose A/ x and A/ 2 are lower and up- 
per boundaries, respectively, of a bandpass containing 
an absorption feature of interest. A "pseudo-" equiva- 
lent width is given by the spectral line index defined, for 
example, by Equation A2 of Cenarr o et al. (2001; see also 
lGonzkledlT99l iCardiel et al.lfl998h : 



V/2 



W(A) = / [1 - S(X)/C(X)]dX, 



(8) 



where C (A) is is an estimate of the continuum flux. Let- 
ting A ri ,A r2 and Af, 1 ,Aft 2 denote boundaries of nearby 
continuum regions redward and bluewar d, respectively, 
of the feature bandpass, C(A) is given by iCenarro et al.1 
(|200lT) : 

C(X) = S b ±-± + S r ^±. (9) 



X r — Xh 



X r — Xh 



Here X r and A& are the central wavelengths of the red and 
blue continuum bandpasses, respectively, and 5), and S r 
are estimates of the mean level in the two continuum 
bands: 

5 b = (A h2 -A bl )- 1 f Xb2 S(X)dX; 

S r = {X r2 -X^)- 1 f ' 2 S(X)dX. 



The Lick syst em of indices (iBurstein et al.l 11984 
iFaber et al.l 119851 : IWorthev etaTI |1994|), whilein wide 
use, is not applicable to our spectra, as the relevant 
Lick magnesium passband alone is larger than our en- 
tire spectral range! Instead we take advantage of the 
high MMFS spectral resolution and define a set of six- 
teen indices, each of which measures the flux in a single, 
resolved spectral line. 

To aid in the determination of feature and continuum 
bandpasses appropriate for the MMFS spectra we pro- 
duced a high-S/N composite spectrum for each galaxy 
by averaging its individual stellar spectra (with relative 
redshifts and probable nonmembers removed). Figures 
[T71 and [TH] display composite dSph spectra for blue and 
red channels. We define feature bandpasses to cover the 
sixteen atomic absorption lines most prominent in these 
spectra. In practice the width of each feature bandpass is 
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Fig. 14. — Comparison of MMFS velocity measurements with previous results. For stars in common with the specified RV sample, 
each panel plots the difference in measured velocity (AV = Vmmfs — Mother) as a function of the MMFS velocity. The bottom subpanel 
within each panel plots internal comparisons from stars with multiple MMFS measurements. Filled squares denote stars for which both 
the MMFS and comparison measurements have quoted uncertainties less than the median uncertainty in their respective sample. Open 
squares/triangles denote stars for which the MMFS measurement has uncertainty less/greater than the median MMFS uncertainty and the 
comparison measurement has uncertainty larger/less than the median from its sample. Crosses denote stars for which both the MMFS and 
published measurements have uncertainties greater than the respective sample median. Histograms in right-hand sub-panels, individually 
normalized for legibility, indicate distributions of AV. 

approximately proportional to the mean line strength in 
the composite spectra. In most cases, blue and red con- 
tinuum bandpasses are adjacent to the feature bandpass 
and span at least 0.25 A (~ 5 pixels). Where possible 
(in the absence of a nearby feature) we use larger contin- 
uum passbands. Feature and continuum passbands for 
the chosen indices appear as shaded regions over the blue 
(red) composite galaxy spectra in Figure [T7] (Figure H5|) . 
Table [6] lists bandpass boundaries for the sixteen MMFS 
indices. 

7.1. Measurement of MMFS Indices 

Further processing of individual stellar spectra is re- 
quired prior to measuring line indices. For dSph and 
globular cluster stellar spectra, we begin with the sky- 
subtracted frames obtained during the data reduction 



steps described in Section 01 We use the DOPCOR task 
and the measured radial velocities to shift wavelength 
solutions such that all spectra have zero redshift in the 
solar rest frame. We then truncate spectra such that all 
cover the range 5141 — 5177 A with linear dispersion 0.071 
A pix -1 (blue) and 0.107 A pix -1 (red). These choices 
of dispersion values minimize pixel re-binning. From the 
resulting spectra, which share a common redshift, range, 
and dispersion, we measure the sixteen MMFS indices 
In order to estimate index mea- 



using Equations IMTOl 
surement errors, we assume the amounts by which the 
measured values W differ from "true" index values W* 
follow Gaussian distributions with variances Oyy . For the 
i th measurement of the j th star, 
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Fig. 15. — Comparison of MMFS velocity measurements with previous res ults, within a window spanning the approximate velocity 
range of dSph member stars. Symbols have the same meanings as in Figure [14] The dotted line in each panel provides the best fit to the 
plotted points. The steepness of the slopes in the Kleyna et al. (2004) and Munoz et al. (2006 w/ MIKE) comparisons are significant at 
the 99.8% and 98.5% levels, respectively. 



TABLE 6 
MMFS Spectral Line Indices 



Index (Jentral Bandpass Kcd Uontmuum Blue Continuum 

(A) (A) (A) 

Fel 5 i42 5141.50 - 5142.00 — 5141.10 - 5141.40 — 5142.05 - 5142.35 

Fel/Nil 5 i4 3 5142.10 - 5143.40 5141.10 - 5141.50 5143.80 - 5144.60 

Nil 5 i 46 5146.00 - 5146.80 5145.65 - 5145.90 5146.90 - 5147.30 

Fel 5 i 48 5147.85 - 5148.45 5147.60 - 5147.80 5148.60 - 5148.90 

Fel 5 i 5 i 5150.40 - 5151.20 5148.90 - 5150.00 5151.25 - 5151.65 

Fel 5 i52 5151.65 - 5152.45 5151.20 - 5151.60 5152.50 - 5152.95 

CoI/TiII 5 i54 5153.85 - 5154.40 5153.55 - 5153.80 5154.55 - 5154.90 

Nil 5 i 56 5155.55 - 5156.00 5155.25 - 5155.50 5156.15 - 5156.40 

Fel 5 i 59 5158.65 - 5159.35 5158.00 - 5158.45 5159.45 - 5159.95 

Fel 5 i62 5161.90 - 5162.70 5161.00 - 5161.60 5162.85 - 5163.50 

Felsiee 5165.90 - 5166.60 5165.60 - 5165.85 5166.65 - 5166.85 

Mgl 5 i67 5166.80 - 5168.00 5164.25 - 5165.25 5169.55 - 5170.70 

FeI/FeII 5 i 69 5168.50 - 5169.45 5168.00 - 5168.40 5169.50 - 5170.00 

Fel™ 5171.30 - 5171.90 5170.70 - 5171.20 5171.95 - 5172.15 

5171.90 - 5173.35 5169.55 - 5170.70 5174.35 - 5175.50 



rci 5 
Mgli 

Til 5 i74" 5173.50 - 5174.00 5173.15 - 5173.40 5174.10 - 5174.50 
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Fig. 16. — For MMFS stars in common with previously published samples, distribution of velocity deviations AV = Vmmfs — V a ther 



normalized by combined errors a = y a \,[MFS + a other- 
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X (A) 



Fig. 17. — Composite spectra produced by averaging (with redshift: 
stellar and twilight spectra acquired with MIKE's blue channel. Black 
line indices. Blue and red shaded regions correspond to the associated 
subset of the sixteen indices. 

As in Equation [2J the values follow the standard nor- 
mal dist ribution, and we ado pt the error approximation 
given by ICardiel et all l|1998ft , 

ci - c 2 W 



aw ~ (3 



(SN) 



(11) 



modified by our introduction of scaling parameter (3. 
Constants c\ and c 2 depend on th e index definit i on an d 
are given by Equations 43-44 of ICardiel et alj (|1998f l. 
and we use our Equation [1] to estimate (SN), the mean 
S/N per pixel. After taking logarithms and replac- 



0.55 (Section 15. ip . we obtain 
2 \ _ i~„ra2\ 



ing log[e? ] with S i: 

log[(W„ - W^f] - log(^) = log(/3 2 ) - 0.55. We es 
timate (3 uniquely for each of the sixteen indices on both 
channels via linear regression, using the approximation 
W*j ~ Yl!i=i WijWij I Y^a=\ w ij (where Wjj is the inverse 
of the variance obtained from Equation QT] with (3=1). 
For the sixteen blue indices, estimates all fall within 
the range 3.45 < (3 < 4.89. For the eight red indices, 



s and spectra from probable nonmembers removed) individual dSph 
shaded regions correspond to feature bandpasses for sixteen MMFS 
continuum bandpasses. To avoid overlap, each panel shows a unique 

2.72 < < 3.95. Among all indices, the resulting errors 
have median values in the range 0.05 < <tw/A < 0.10. 

7.2. Individual Index Results 

The upper-left panels in Figures fl9l and [20l display, for 
each index, the measured value as a function of the mean 
S/N per pixel (calculated over the range 5141 — 5177 
A). Because the red channel has poor sensitivity toward 
the blue end of our spectral coverage (see the composite 
red spectra in Figure IT8|) , we consider only the eight red 
indices with feature bandpass redward of 5158 A. For a 
given index, the width of the distribution of W values 
correlates with both (S/N) and the width of the feature 
bandpass 

The remaining panels in Figures flUl and [2"01 display em- 
pirical relationships between individual MMFS indices 
and other stellar parameters: V-I color, radial velocity, 
and Mgl5i73 index. For several indices there is a clear 
trend toward higher index values at redder colors (upper- 
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right panel); this behavior indicates the expected de- 
pendence of line strength on effective temperature. The 
lower-left panels plot W against the measured radial ve- 
locity; for clarity in these panels the plotted points corre- 
spond exclusively to spectra from Carina targets. Along 
the radial velocity axis, most points are clumped in rel- 
atively narrow distributions about Carina's systemic ve- 
locity (~ 225 km s _1 ). That the velocity outliers with 
respect to the narrow Carina distribution tend also to 
be Mgl index outliers indicates the ability of the Mgl in- 
dices to distinguish dSph members from interlopers. We 
expect interlopers with large Mgl line strengths to be 
due to Milky Way dwarf stars with high surface gravity. 
Thus the spectral indices provide a basis for dwarf/giant 
discrimination similar to that derived from photometric 
techniques that use narro w-band filters (e.g., DDQ51) 
sensi tive to surface gravity ()Geislerl ll984: Maiewsk i et alJ 

Hoop). 

Finally, when individual indices are plotted against 
one-another, indices corresponding to the same ele- 
ment (s) display monotonic, linear relationships with 



slopes approximately equal to the ratio of feature band- 
pass widths (if indices are stated as a fractions of their 
bandpass widths, these relations tend toward 1:1). For 
example, the lower-right panels of Figures [TO1 and [2U1 plot 
each index against the Mglsi73 index. In this case the 
Mgi5i67 index exhibits the expected monotonic, linear 
relation. Noting the behavior of the Fel indices, it is 
immediately striking that many exhibit a bifurcation. 
With respect to the trend followed by stars with small 
and intermediate Mglsi73 index values, stars with the 
largest Mgl5i73 values have systematically small Fel in- 
dices. Such stars tend to be the outliers in the radial ve- 
locity and Mgl index joint distributions (lower-left pan- 
els); thus the bifurcation is due primarily to the presence 
of interloping stars, and dSph members generally follow 
a linear relationship in the iron-magnesium index plane. 

7.3. Composite Indices 

We have noted that indices for different absorption 
lines corresponding to the same element (s) correlate lin- 
early. We can improve precision by computing their sum. 
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Fig. 19. — Measured index value vs. spectral and stellar parameters for blue-channel spectra. Each panel corresponds to one of the 
sixteen MMFS indices. Points in the lower- left panel correspond to spectra from only Carina target stars (interlopers included), and show 
a clear separation in Mgl indices between dSph and interloper populations. 



For blue-channel spectra we define a composite Fe index: 
EFe blue = 0.160(FeI 5 i42 + Fel/Nil 5143 + Fel 5148 (12) 
+FeI 5 i5i + Fel 5 i 52 + Fel 5 i 59 + Fel 5 i 62 
+FeI 5 i 66 + FeI/FeII 5 i69 + FeFjm), 
where the sum is scaled by the ratio of the Felsi43 band- 
pass width to the sum of all Fel bandpass widths (this 
scaling effectively weights each index by its bandpass 
width). Similarly, for red-channel spectra the compos- 
ite Fe index is given by 

EFe rcd = 0.253(FeI 5159 + Fel 5 i 6 2 + Fel 5 i 66 (13) 
+FcI/FeII 5 i69 +FeI 5 m), 
where in this case the straight sum is scaled by the ratio 
of the Felsieg bandpass width to the sum of the included 
Fel bandpass widths. For both blue- and red-channel 
spectra we compute the composite Mg index as 

EMg = 0.547(MgI 5167 + Mgl 5173 ), (14) 

where the scale factor is the ratio of the Mgl.5173 band- 
pass width to the sum of the two Mgl bandpass widths. 



We compute errors a-sw associated with composite in- 
dices by adding in quadrature the estimated errors aw 
of the contributing individual indices, and then scaling 
by the constant multipliers in Equations [T3TfT4l Figures 
[2T1 and [22] plot for composite indices the same relations 
shown for individual indices in Figures [19] and [20] The 
composite indices show the same features and behaviors 
as the relevant individual indices, only with less scatter. 
Figure [23] plots estimated errors as a function of mean 
S/N per pixel. 

We note that EFeblue includes a small contribution 
from a Nil absorption line at 5143 A. This line is blended 
with two Fel lines in the MMFS spectra. Because the 
Fel/Nil5i43 index behaves similarly to the other Fel in- 
dices (see Figures [T5] and [2U]) and because it has among 
the highest S/N of any individual Fe index, we choose 
to include it in EFcbiuc- For similar reasons we choose 
to include the FeI/FcII 5 i 6 9 in both EFeblue an d SFe ro d 
despite a small contribution from a blended Fell feature. 
Exclusion of either index from the composite Fe index 
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Fig. 20. — Same as Figure IT9l but for red-channel spectra. 

does not significantly alter any of the correlations or fea- 
tures seen in Figures [21] and [22l nor the distribution of 
errors seen in Figure 1231 

For stars with multiple index measurements, Figure l24l 
plots the distribution of AEVF = EWjy — "SWj, normal- 
ized by the associated errors, where TiWj is the weighted 
mean. The best-fitting Gaussian distributions (solid 
curves) are nearly identical to the standard normal dis- 
tribution (dotted curves), indicating that the estimated 
errors are valid. 

8. CONCLUSION 

We have described the instrumentation, data acqui- 
sition and data reduction procedures for a new spec- 
troscopic survey that targets large numbers of stars in 
nearby dSphs. In forthcoming papers (Walker et al. in 
preparation) we present the entire data set described 
herein and provide kinematic analyses that consider both 
equilibrium and tidal interaction models. We use the 
spectral indices measured from stars in seven globular 



clusters to calibrate relations between the MMFS spec- 
tral indices and iron abundance, and to identify dSph 
interlopers independently of velocity. In the end we have 
velocity and line-strength measurements for nearly 3800 
probable members in the four observed dSphs. Alone 
and/or in combination with data published from other 
contemporary spectroscopic surveys, the MMFS data set 
enables us to map the two-dimensional behavior of dSph 
stellar kinematics and chemistry with unprecedented pre- 
cision. 
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